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Abstract 

Correlation functions in the restricted primitive model are calculated within a field-theoretic 
approach in the one-loop self-consistent Hartree approximation. The correlation functions exhibit 
damped oscillatory behavior as found before in the Gaussian approximation [Ciach at. al., J. 
Chem. Phys. 118, 3702 (2003)]. The fluctuation contribution leads to a renormalization of both 
the amplitude and the decay length of the correlation functions. The renormalized quantities show 
qualitatively different behavior than their mean-field (MF) counterparts. While the amplitude and 
the decay length both diverge in MF when the A-line is approached, the renormalized quantities 
remain of order of unity in the same dimensionless units down to the coexistence with the ionic 
crystal. Along the line of the phase transition the decay length and the period of oscillations are 
independent of density, and their values in units of the diameter of the ions are 1 ~ 1 and 
2n/ai « 2.8 respectively. 



* Dedicated to Prof. R. Evans on the occasion of his 60th birthday 



I. INTRODUCTION 



Distribution of charges in ionic systems has been a subject of considerable interest for 
many years. In theoretical studies of systems such as molten salts, electrolytes or ionic liquids 
the interaction potential is often approximated by the restricted primitive model (RPM), 
where identical hard spheres carry positive or negative charges with equal magnitude. Earlier 
studies concentrated on correlation functions for charges in ion-dilute systems 111- 2J, y], and 
later the ion-dense systems were also investigated However, the 

question of the form of the correlation functions in the liquid (fused salt) at the coexistence 
with the ionic crystal remain open. Whereas the correlation functions in a liquid can be 
quite accurately described within liquid theories such as the mean-spherical approximation 
(MSA) , the transition to a crystalline phase is not predicted by the liquid theories. 
The liquid-solid transition in the RPM was determined within the density functional theory 
of freezing in Ref. [10], but the charge-charge correlation function is not discussed in this 
work. Both the phase transition and the correlation functions can be found on the same 
level of approximation within the framework of the field theory for the RPM introduced in 



Ref. |l lj| . In this work we calculate the correlation functions in the RPM by applying the 
Brazowskii-type approximation. We pay particular attention to the decay length, the period 
of damped oscillations and the amplitude along the phase coexistence with an ionic crystal 
that was recently found in Ref. |l2j |. 

The behavior of the charge-charge correlation function hr>(x) depends on thermodynamic 
conditions. In Fig.l a portion of the density-temperature phase diagram, as obtained in 



different theoretical 2L 3f and simulation 



13| studies, is shown schematically. The solid 



lines represent boundaries of stability of the uniform fluid phase, and regions corresponding 
to different qualitative properties of ho(x) are separated by the dashed lines. 

The Kirkwood line (K) separates the high-temperature/low-density region corresponding 
to an asymptotic monotonic decay of correlations, from the low-temperature/high-density 
region corresponding to exponentially damped oscillatory behavior of ho(x). In theories of 
the mean-field (MF) type the Kirkwood line is found to be a straight line S = Sk = const, 
where 

s = -. (1) 

Pi 

In the above p$ is the dimensionless number-density (volume measured in <r 3 units, where a 



2 



p 



FIG. 1: A portion of the density-temperature phase diagram for the RPM as obtained in MF, 
shown schematically. Solid lies separate the uniform fluid from the two-phase regions. Dashed 
lines denoted by 'K' and 'A' represent the Kirkwood line and the A-line respectively. T and p are 
in arbitrary units. 

is the ion-diameter), and 

* DakT 

T = =— (2) 

e z 

is the temperature in standard reduced units (D is the dielectric constant of the solvent and 
e is the charge). In different theories the value of Sk is somewhat different, and typically 
Sk ~ 10. There is no controversy concerning this line, but there are no general reasons for 
this line to be straight beyond MF [3J. 

When the A-line is a ppr oached from the high-temperature side, the decay vent of corre- 
lations diverges HQ 15]. The amplitude of the correlation function diverges [16] as well 



on the MF level of the field theory |ll(. The low-temperature side of this line corresponds 
to a charge-ordered structure, where the nearest-neighbors (nn) are oppositely charged, and 
the A-line represents a continuous transition between the uniform and the charge-ordered 
phases. The slope of this line shows a strong dependence on approximations and assumptions 



made in different theories 



13. In particular, in the MF approximation for the field theory 



the A-line is S = S\ ~ 1.6 [TJ|, whereas the MSA yields S\ = 0. The charge-ordered 
liquid phase was observed neither experimentally nor in simulations. A natural conclusion 
would be to question validity of theories that predict the existence of the A-line. However, 
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ihe A-line of continuous transitions to the charge-ordered phase was observed in simulations 
l?! Q, of the lattice version of the RPM (LRPM), in agreement with predictions of 
the lattice version of the field theory Some other theories QHQ yield for the 



LRPM similar results. The MSA, however, predicts an absence of such a transition in the 
LRPM. One might argue that the MSA theories are correct only in the continuum space, 
and the field theory yields correct results only on the lattice. This is indeed the case for the 
field theory in the MF approximation. Beyond MF, with the charge fluctuations accounted 



for in a Brazovskii-type approximation \l'A |24j, |25j|, the decay length of the charge-charge 
correlation function is found to be finite in the continuum-space RPM Q|. Thus, when 
the fluctuations are included in the field theory, there is no line of continuous transitions, 
in agreement with the MSA. At the same time the first-order transition to the ionic crystal 
is found Q for the range of densities and temperatures that agree with simulations Q. 
Thus, within the field theory we can calculate the correlation functions at the coexistence 
with the ionic crystal. 
The field-theory results 



24, also shed light on the nature of the A-line that appears 



in MF. Beyond the Kirkwood line there is no qualitative change of the average structure of 
the liquid phase, which is represented by the correlation function. On the other hand, the line 
S ~ 1.6 separates the phase diagram (/?q,T*) into two regions that correspond to different 
forms of the most probable instantaneous distributions of ions. At the high-temperature 
side of the line S = S\ the distribution of ions is most probably random. For S < S\, 
however, the most probable instantaneous distributions of the ions correspond to charge- 
density waves with a wavelength ~ 2.5a, where a is the hard-sphere diameter. According to 
the above prediction, in real-space representation charge-ordered clusters with a nn distance 
~ 1.27a should dominate over randomly distributed ions for S < S\ in instantaneous charge 
distributions (i.e. in microscopic states). Such structures were indeed seen in simulation 

nn 

snapshots |26J, |2j| . In MF the average density distributions are approximated by the most 
probable distributions. When in a particular system the average- and the most-probable 
distributions agree on a qualitative level, the MF results are qualitatively correct. This is 
not the case in the RPM, therefore the line S = S\ is erroneously interpreted as a line 
of continuous phase transitions in MF theories. On the other hand, the occurrence of the 
A-line on the MF level is an indication of the first-order transition to the ionic crystal, 
because beyond MF the first-order transition is found in the phase-space region where the 
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charge-ordered individual microstates are more probable than the microstates with randomly 
distributed ions. 

In this work we focus mainly on the charge-charge correlations for S < S\, and calculate 
the number-number correlation function as well. In Ref.|16| the correlation functions are 

n 

calculated within the field-theoretic description developed in Ref. |11|] only in the Gaussian 
approximation. As already mentioned, at the A-line the decay length and the amplitude of 
the charge-charge correlation function both diverge. Our purpose is to verify if the instanta- 
neous order has some effect on the correlation functions beyond the Gaussian approximation. 
The results of Ref. are extended beyond the Gaussian approximation within the Bra- 



zovs 
Ref. 



rii theory 



281 ] . It turns out that at the coexistence with the ionic crystal, obtained in 



12 1 , the amplitude and the decay length are quite small. Our results show that large 



decay length and amplitude of the charge-charge correlation function can be found only in 
the metastable liquid. 

II. FORMALISM 



A. Field-theoretic description of the RPM 

In the field theory derived for the RPM in Ref. [ill. 



25 



29j | the fluctuating fields are the 
local deviations from the most probable number- and charge densities, r/(x) = p*(x) — p£ = 
P+( x ) + pl(x) — Pq and 0(x) = p+(x) — pl(x) respectively. The most probable number 
density of ions is denoted by Pq, and the subscripts + and — refer to cations and anions 
respectively. Asterisks indicate that all densities are dimensionless (the unit volume is a 3 ), 
(ft is the charge density in e/cr 3 units. We focus on systems that are globally charge-neutral, 



fx) = 0, 



(3) 



where in this paper we use the notation / = f c?x. The fields 77 (x) and 0(x) are thermally 
excited with the probability density [ill GI 1^0 1 



p[<f>,rj\ 



exp(-/3Afi MF [0,r/]), 



where 



Dr) / D(f>e 



(4) 



(5) 



is the normalization constant and Afl MF [4>,r]} = Q MF [4>, p*} — Q mf [0,Pq], where for 6 = 
and p* = Pq the Q MF [(f), p*} assumes a minimum. In the theory studied in Refs. 
Q MF is approximated by 

Q MF [<P,p*} = F h [<p,p*] + U[<j>] -ill p(x). (6) 

J X 

In the above p is the chemical potential of the ions, Fh[4>, p*] = f /ft(0(x), p*(x)) is the 
hard-core reference-system Helmholtz free-energy of the mixture in which the core-diameter 
a of both components is the same. We use the Carnahan-Starling (CS) form of fh(4>, p*) in 
the local-density approximation, as in Ref.Q). Finally, the energy in the RPM is given by 



W] = y I I 9(l 1x-l\ - ^i*') = t I Wvtm-K), (7) 

•J x J II J k 



where J k = J <ik/(27r) 3 , and /(k) denotes a Fourier transform of /(x). Contributions to the 
electrostatic energy coming from overlapping cores are not included in (J7|). Here and below 
the distance is measured in a units, and (3* = 1/T* with T* defined in Eq.(J2J). The Fourier 
transform of V(x) = 6(x — l)/x, where x = |x|, is 

V(k) = Air cos k/k 2 , (8) 

where k is in a^ 1 units. 

In the field theory the charge- density fields are considered. The correlation function for 
the fields at points x and x', called the charge- density correlation function, is given by 

G*(|x - x'|) = <0(x)0(x')> = E- 1 Jd V J D0e-^ MF 0(x)0(x'). (9) 

The charge-density correlation function defined in Eq.(J2} is related to the analog of the total 
charge-charge correlation function Hd{x) H ||3|, and in Fourier representation this relation 
is given by 

^M = h D (k) P * + i. (io) 



B. Weighted-field approximation 

Because AQ MF [0,r]} is strictly local, the field 77 can be integrated out exactly |3l| and 
Eq.Q takes the form 

GW(|x - x'|) = H- 1 J D0e-^//M0(x)0(x'), (11) 



where 

m^W\ = U f 0(x)C( x - x >(*') + E / ^ m w- (12) 

In the above 

C° (x, x') = C^(x, x') + / 8 *(x - x'), (13) 
where in Fourier representation C^L is given by 

CUk) = -J^i . (14) 

The fa and the coupling constants Aim are functions of pj$, whose exact forms depend on 



fh{4>,P*) j3J|, and have not been determined yet. Approximate expressions for Aim can be 
obtained in the weighted field approximation (WF) introduced in Ref . Il ll . and described 

nnn 

in more detail in Refs. |25l . I29L l30j |. In the WF the field r/(x) is approximated by its most 
probable form for each given field 0(x). Another words, for a given field 0(x), the field i](x.) 
is determined by the minimum of /3AQ MF [<p, 77], i.e. it assumes the form 



W f(0(x)) = ^^(x) 2 ", (15) 

n 

where a n are functions of Pq Q; 13 ■ The average density of ions in the WF is given by 

(p*(x)) = pi + J L>0e-^//M WF (0(x)). (16) 

The WF approximation is equivalent to the lowest-order, zero-loop approximation for the 
exact TC e ff[4>] jsi|- In the WF fa = j^lj], and ^4 ?T1 and a„, given in terms of derivatives of 
fh{<t>, p*) at (f) = and p* = p* Q , are functions of Pq |2S|- 

In this work we limit ourselves to the 6 theory, with TC e ff approximated by 

FHwfW = \j I *(x)C$,(x - x')0(x') + ^ / 4 (x) + ^ /* 6 (x). (17) 

*/ x */ x' */ X */ X 

The explicit forms of A^pj^O and Ao(po) in the WF theory with the CS approximation for fa 



are given in Appendix A [12j. In the density interval 0.9 < Pq < 0.1541 the above functional 
has to be supplemented with a term oc <f) 8 , to ensure stability of PT-Cwf^] for <fi —>■ oo. In 
this work we consider only pj > 0.1541; for this range of densities Aq{Pq) > 0. One should 
remember that the accuracy of the results obtained with the above functional depends on 
how large are the neglected contributions to /3Ti, e ff[4>]. 
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In the 6 theory the average density in the uniform phase, Eq.([16jh is approximated 



by 



(p*(x))=p* + a 1 (0(x) 2 ) + |(0(x) 4 ), (IS) 



where the explicit forms of ai(po) an d °2(Po) ? obtained in Ref. [12], are given in Appendix 
A. On the same level of approximation the number- density correlation function for x^x' 
is given by 



4 

where 



C w (x,x') = ( WF (x) WF (x')) co " = a?(0 2 (x)0 2 (x')) co " 
a\a 2 ' 



(V (x)0 2 (x')> co " + (x <- x')) + ^(0 4 (x)0 4 (x')> COT \ (19) 



(/(0(x))^(0(x'))) co " = </(0(x)M0(x'))> - (/(0(x)))(^(0(x'))). (20) 



C. One-loop self-consistent Hartree approximation 



Even for the approximate effective Hamiltonian (1171), given in Eq. (jllj) cannot be 
calculated analytically. In the perturbation theory G^ is given by Feynman dia- 



grams 



331 ] with the 2n-point vertices Azn- The vertices at x and x' are connected by 



lines representing G^(x — x'), where G° is inverse to C° defined in Eq.()14|). and all lines are 
paired. The corresponding expressions are inte grat ed over all vertex points, or in Fourier 
representation over all G?|L (fc)-line loops, where jlll ll^ 



G° H> 1 (k) = Cl(k)=(3*(S + V(k)). 



Note that 



where 



S 



47T 

x D 



x D = k d g = a/ 47r/3*Po 

is the dimensionless inverse Debye length. 

The one- loop contribution to (Fig. 2a) is proportional to A^Qo, where 



(21) 



(22) 



(23) 



So 



GUk)- 



(24) 



S 



^ 



a be 

FIG. 2: Feynman diagrams contributing to to two-loop order. The shaded circle and the 
bullet represent A4 and Aq respectively. Lines represent G^,. In the self-consistent theory the 
lines represent G^. 

Note that Eq. (j21j) can be written in the form 

C° H (k) = (3*{r + AV(k)) (25) 

where 

P*r = \ + (3*V(k b ), (26) 
Po 

k = kf, corresponds to a minimum of V(k), and 

AV(k) = V(k) - V(k b ) ~ k ^ kb v 2 (k - k b f + 0((k - k b ) 3 ). (27) 

The integrand in Eq. (j24j) diverges for k = k\, when (3*Tq = 0. For k — > kb Eq. (|24jl is of a 
similar form as the corresponding integral in the Brazovskii theory [28]. and Go diverges for 

ler contribution to C^ik) is given 



(3*Tq — > 0. In the effectively one-loop 6 theory ()17j) . anot 
by a diagram (Fig. 2c) that is proportional to AnGn 25l. 



32 



33 1 . 



29(. The symmetry factors of the 



graphs are calculated according to standard rules 

In this work is calculated in the effectively one-loop approximation, i.e. as in Ref. 
the diagrams shown in Fig. 2 a,c are taken into account. In the self-consistent, effectively one- 
loop Hartree approximation the inverse correlation function assumes the form 



where 



and 



Ci(k) =r + f3*AV(k) (28) 



r = CM=?*r + ^ + ^pl, (29) 



^(r o ) = (0(x) 2 ) = ^G^). (30) 






b, 



b 2 



FIG. 3: Feynman diagrams contributing to G^(xi,X2) given in Eq. (|19j) . Open circles represent 
external points xi and X2. Lines connecting points x and x' represent G^(x — x'). The black box 
represents 1/70,2 , the vertex from which the box and the two thin lines emanate represents —72,1 
and the circle with a cross inside represents 1 (see Appendix A, Eas. (|47|) and (|51|)V Symmetry 



factors are calculated according to standard rules 



32, 



33 



The integral (|30j) diverges because of the integrand behavior for k —>■ 00. However, the 
contribution from k —>■ 00 is unphysical (overlapping hard cores). When the fluctuations 
with k ~ kb dominate, which is the case for r -C P*v 2 kl, then the main physical contribution 
to Q(ro) comes from k « kj,- In this case 



C»(k)^r + (3*v 2 (k-k b f 



(31) 



and the regularized integral is [2i 

1 



k r + (3*AV(k) 



-ro— >0 



2aVT 



r + (3*v 2 (k- k b y 



(32) 



where 

a = k 2 b /(An^ 2 ). (33) 

The self-consistent solution of the pair of equations (|29|) and ()32|) can be found analyti- 
cally. The explicit expression for ro is given in Appendix B. 

The above results allow for obtaining the number- density correlation function G m at the 
same level of approximation. From Eq. ()19j1 we obtain 

G^(x,x') = {2a\ + \2a x a 2 Q{r G ) + 18a 2 £(r ) 2 )G%(x, x') 2 + 3!a 2 G%(x, x') 4 , (34) 

with the relevant dia gram s shown in Fig.3, and the symmetry factors calculated according 



to standard rules 



32, 



33). 
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III. PROPERTIES OF THE CORRELATION FUNCTIONS 



A. Gaussian approximation 

In Ref . [r| the correlation function G^ was studied on the level of the Landau- type the- 
ory, where is approximated by G^. In Fourier representation the above approximation 
for the charge- density correlation function is given in Eq.(|21|l{llHl3|. The Fourier transform 
G^{k) of the correlation function defined in Eq. ffTTj) reduces to G^Jk) given in Eq. (j2*T|) when 
all ^-dependent terms in j3AQ MF , except from the Gaussian part | f k 4>(k)C^(k)(f)(— k), are 
neglected. 

Let us first focus on the behavior of G°,Ak) for k — * 0. From Eqs. (|2~Tj) and (JHJ) we obtain 



GUk) fc 2 



Po k 2 + x 



2 ' 

D 



(35) 



where the dimensionless Debye length xd is given in Eq. (j23|) . Note that Eq. ()35|) (with 
5q. (|Tn|) ) is consistent with the exact second- moment condition of Stillinger and Lovett 
and with the Debye-Hiickel result. 
For finite values of k the properties of G^(k) agree with the results of other theories 2 
only qualitatively The qualitative agreement with the GMSA results for hoik) 5| is 
found only in the phase-space region given by 

S^S x = -V(k b ), (36) 

where 5* is defined in Eq.([T]). The straight line S = S\ is equivalent to r = 0, and represents 
the boundary of stability of the functional Tiy/F m the (pq, T*) phase diagram jl2, 25]. At this 
line the second functional-derivative of H.wf vanishes for the most probable fluctuations 25, 

. The line of instability of riwFyP] is of course the same as the line of instability of the 
functional AQ MF (see Eq.flHJ) and Refs.QHQ)- As shown in Ref. Q, the decay length 
and the amplitude of G^^x) both diverge for S — > S\. 

The charge- density correlation function in real-space representation is obtained in Ref. 
by a pole analysis of G°^{k). The Kirkwood line in the Landau theory [3] is 5* = Sk, where 

47rcoshatf 2 

Sk = n , tanha^- = — , (37) 

a K a K 
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HQ 



and Sk ~ 11.8. In the GMSA and in some other MF theories [5|, [15( the Kirkwood line is 
also a straight line, but in the GMSA its slope was found to be Sk = 47r/x^ ~ 8.33 [5|. For 
S > Sk there are two inverse decay- lengths, a\ and a 2 , and 

xG^(x)=A^e^' + A^e^". (38) 

For S — > S K the two imaginary poles of G^Jk), ia\ and ia 2 , merge together and for S < Sk a 
pair of conjugate complex poles, zcto±ai, appears. As a result, G^(x) exhibits an oscillatory 
behavior for S < Sk, 

xG%(x) = -A H sin( ai x + 9)e~ a ° x , (39) 

with the amplitude A^, the phase 9 and the inverse lengths a and «i given in Ref.[3|. 
When the A-line is approached, both a and A^ diverge as oc (S — Sx)^ 1 ^ 2 , whereas 9 oc 
{S - 5 A ) 1/2 and a x -> k b . 

The number-density correlation function in the above Landau theory for the RPM has 
the trivial form 

<4 = 7o>Kx-x'), (40) 

because the microscopic structure associated with hard spheres in the coarse-grained de- 
scription is suppressed. 



B. Charge-density correlation function in the self-consistent one- loop approxima- 
tion 

In this section we consider the charge-density correlation function Q in the approxima- 
tion outlined in sec. 2. 3. As already explained, our Brazovskii-type approximation relies on 
two important assumptions, and therefore it is valid only in a limited region of the phase 
diagram. In this section we determine the explicit form of the correlation functions and 
also the region of validity of our results. The approximate expression ()28|) for the inverse 
charge- density correlation function can be rewritten in the form 

c H{ k)=(s» + ^y. (4i) 

This formula is analogous to Eq. (j21j) in Ref.Q], except that S = T* / p* Q is replaced by 

S H = T*r -V(h). (42) 
12 



l_ 

10 

* 

T r 



15 



20 







FIG. 4: The inverse decay length ao of the charge-density correlation function, as a function of 
T*ro(/3g, T*), where ro(po,T*) is the self-consistent solution of Eas.(|29|) and For the explicit 

form of ro see Appendix B. The line cko^* 7 ^) splits into two branches representing a\ and a2 at 
the Kirkwood line T*tq = Sk + V{kb). Dashed line is the inverse decay-length obtained from the 
approximate form of C^{k) (Ea.(|31|)). 



The charge- density correlation function in real space can be obtained by following the 
pole analysis of Ref.|l6j|. just by substituting S H for S. The inverse lengths a and oc\ are 
shown in Figs. 4 and 5, where the numerical results obtained in Ref.fl3| are used. Note that 
at the A-line corresponding to the infinite range of correlations (ao = 0), we find T*r^ = 0. 
From the explicit result for r (Appendix B) it follows that r > for T* > 0. For T* > 
the range of charge- density correlations is thus finite, in agreement with the results of the 
MSA and related approximations 

as 

, as well as with simulations. 

The dashed line in Fig. 4 represents ao for the approximate form of C^(fc) given in 
Eq.()31j). The approximate form of Q, Eq. (j32|) . and in turn the form of r are based on the 
above approximation for C^k). As seen form Fig. 4, Eq. (j31j) yields a good approximation 
for ao when T*r < 1 jl^ . 

Note that the approximation for Q(tq) is valid only for r <C fi*Vik\- For r 3> /3*V2kl 
the fluctuations with wavenumbers other than k ~ fc& give a relevant contribution to the 
integral in Eq. (j3H|) . and its value depends significantly on the regularization procedure. In 
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T r 



o 



FIG. 5: ot\ = 2n/\, where A is the period of damped oscillations of function of 

T*ro(pQ, T*), where tq(pq,T*) is the self-consistent solution of Eas. (|2§|) and and its explicit 
form is given in Appendix B. 

addition, in the Brazovskii 4 -theory the two-loop (Fig. 2b) and higher-order diagrams can 
be neglected if A±y/ (3*V2 <C r \2s^ . From the two above conditions we obtain the region in 
the phase diagram, 

v 2 kl 



T* < 



A^ 



(43) 



where Eqs. (129)1 and (|3*2"|) give accurate results in the field theory corresponding to the Hamil- 
tonian ()17j) . The line T* = V2k^/A\ is shown in Fig. 6 together with the coexistence line 
obtained in Ref . Q| > and the K-line (j44j) discussed below. Note that the assumptions of our 
theory are satisfied near the phase coexistence for the high densities pg > 0.5. 
In the approximation derived in sec. 3 the Kirkwood line is given by 



S H (p*,T*) =S K *tll.8, 



(44) 



where S is given in Eq . (J42j) . and the explicit expression for Tq(pq, T*) is given in Appendix 
B. The line f!44j) is represented by the dotted line in Fig. 6. Note that this line lies very far 
from the region of validity of the Brazovskii approximation. Note also that from Eq. (j31j) 
one obtains a\ = The approximation (|31|) fails of course when the Kirkwood line is 
approached, because at the Kirkwood line and beyond a,\ = 0. We conclude that the 

14 



0.2 0.3 0.4 0.5 , 0.6 0.7 0.8 




FIG. 6: Thin solid line represents T* = v 2 kf/Al- For T* < v 2 k\IA\ the results of the self- 
consistent one- loop approximation with ()3U|) approximated by Q32JI are correct. In the rest of the 
phase diagram either the remaining Feynman diagrams (including the one shown in Fig. 2b) have 
to be included, or the approximate form of G(tq) (|32|) is not valid. Thick solid line represents 
the phase transition found in Ref.^J. Dotted line is the Kirkwood line (|44|) in the self-consistent 
one-loop approximation. 

assumptions of our version of the Brazovskii theory are satisfied near the phase coexistence 
for Pq > 0.5. When the K-line is approached, the theory becomes less accurate, and it fails 
completely beyond the K-line. On the other hand, beyond the K-line the simple Gaussian 
approximation yields results consistent with the exact predictions (see (|55|)). 

C. Correlation functions at the coexistence with the ionic crystal 



The transition between the liquid and ionic crystal was found in Ref. 12j in an approach 
consistent with the self-consistent one-loop approximation outlined in sec.2C. Note that 

15 



0.3 0.4 0.5 0.6 0.7 0.8 

Po 



FIG. 7: ^o(pq) along the coexistence line T*(pq) between liquid and the charge-ordered phase in 
the (/> 6 -theory. 

in the one-loop approximation T*C^{k) depends on the thermodynamic state through a 
single variable S H . Since depends on the thermodynamic state only through S H = 

T*r + S\ (Eq.()42|)). we need to know T*r at the coexistence. From the explicit result 
for ro (Appendix B) and the result for the coexistence line T*(pq) obtained in Ref . jisjl . we 
obtain vq(pq) shown in Fig. 7, and in turn we find T*r m 0.2 for 0.25 < Pq < 0.8. The 
inverse lengths along the coexistence are thus nearly constant and independent of density, 
and their numerical values are found to be a® ~ 0.97 and oil ~ 2.26 on the basis of the 
results obtained in Ref. [l^ . The approximate results for the amplitude and the phase are 
Q -A H « 0.23S H /VS B - S\ « 0.96 and 9 w 0.88V^ - S x w 3.81 respectively The 
correlation functions along the phase coexistence are shown in Figs. 8 and 9. 

The inverse lengths «o and a,\ depend on the thermodynamic state through a single 
variable also in the MSA and related approximations 5], and xd = kd<j = \J An / S is 
usually chosen. Note that in our theory S H reduces to S in MF; thus, the parameter 



a/4-7t / 'S H reduces to Xjj. At the transition An / S H ~ 2.25. The characteristic inverse 
lengths in the GMSA for xd ~ 2.25 can be read off from Fig. 4 of Ref. [5], and they are 
a ~ 1.5 and ol\ ~ 2.1. More precise values are given in Ref.0 for = 4, where a Q = 1.32 
and aci = 2.89, and for xo = 1.5, where «o = 2.15 and a,\ = 1.22. Rather good agreement 
between the two theories is obtained, although in our theory the decay length is larger than 
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FIG. 8: Charge-density correlation function j3*G^{x) in the liquid phase for S « 1.8 corre- 
sponding to the coexistence with the charge-ordered phase. /3*G^(x) is dimensionless and x is in 
a units. 

in the GMSA. Similar fair agreement is obtained for the amplitude and the phase-shift. 

Note that the form of P*G^ (j) (x) is the same along the phase coexistence, for a large 
range of density. This result is consistent with the structure of the ionic crystal coexisting 
with the liquid, where the nn distance is also independent of density on the same level 
of approximation. This result reflects the fact that the tendency for ordering in periodic 
structure follows directly form the form of the Coulomb potential. 

The number-density correlation function is shown in Fig.9. In our theory the contri- 
bution to the correlation functions associated with hard-sphere ordering is neglected. The 
number-density correlation function shown in Fig.9 results from the coupling between the 
number- and the charge- density fluctuations, and is induced by the charge-ordering, which 
in turn follows from the form of the interaction potential. The exact form of the number- 
density correlation function should contain the entropic contribution due to the hard-spheres 
ordering, in addition to the interaction-induced contribution shown in Fig.9. 



17 



0.05 



0.04 - 




FIG. 9: Number-density correlation function /3*G^ r) (x) (Eg. 1)34 jl ) in the liquid phase at coexistence 
with the charge-ordered phase for p$ = 0.7. Note that the number-density correlation function is 
three orders of magnitude smaller than the charge-density correlation function. 

IV. SUMMARY 

In this work the fluctuation contribution to the charge- and number-density correlation 
functions in the RPM is calculated within the Brazovskii-type theory. In MF the am- 
plitude and the correlation length both diverge when the A-line is approached. At the 
low-temperature side of the A-line an instantaneous order occurs in the majority of the indi- 
vidual microscopic states, i.a. clusters with oppositely charged nn dominate over randomly 
distributed ions. Averaging over microscopic states with somewhat different order leads to 
a homogeneization of the structure downhomogenizationtion-induced first-order transition 
to the ionic crystal. Our purpose here was a verification if the instantaneous order found 
on the low-T side of the A-line leaves some traces in the correlation functions. Our results 
show that the correlation length and the amplitude of G^x), both divergent in MF at 
the A-line, show no indication of the MF singularity when the fluctuation contribution is 
included. Both quantities are of order of unity in the whole stability region of the fluid 
phase. 
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V. APPENDICES 

A. Coefficients A^,Aq, a \ and 02 in the WF approximation 



. 9 (~72,l) 

A4, = 74,0 - 3 , 

70,2 

a _„ 1K (-72,l) (-74,l) 1K (-72,l) 3 (-70,3) /1K (-72,2)(-72,l) 2 
^6 — 76,0 — 10 10 o 40 2 : 

70,2 70,2 70,2 

where 

d 2m + n f h {<p,p) 

In the CS approximation they assume the explicit forms 

1 - 20s + 10s 2 - 4s 3 + s 4 



Aa 



p^(l + 4s + 4s 2 -4s 3 + s 4 



and 



p£ 5 (l + 4s + 4s 2 -4s 3 + s 4 ) 5 ' 



where 



W{s) = 3 - 84s + 360s 2 + 2644s 3 + 1701s 4 - 8736s 5 
+11240s 6 - 8304s 7 + 3861s 8 - 1164s 9 + 240s 10 - 36s 11 + 3s 12 . 



The coefficients a\ and 02 in Eq. (jl8j) are 

72,1 a 2 



r° 

2170,2' 2 4!7o i2 
where 



rji = 74,1 



6(-72,2)(-72,i) 3(-7 3i o)(-72,i) i 



2 



70,2 75,2 
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B. Explicit expression for ro 



By solving Eqs. ()29|) and (|32|) we obtain in the case of To < the explicit expression for 



1 r 



2 



y/u + v - p/3 + y/2(h-(u + v)/2-p/3) + r 
where the following notations are introduced: 



(52) 



u 

h 
h 
P 
Q 



-I(g + 8^ ), , 2 = (_|) 3 + ^r^ -| 

[((« + u)/2 + p/3) 2 + 3{u- vf /A] 1/2 
- (r 2 + 2a 2 T*A) , 



a A T* 2 Ai 



t and a are defined in Eqs.(|2bjl and respectively. 
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